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S 1 Abstract 

q We derive the expressions needed to interpret experiments relating to in- 

O ' 

terplane magnetic coupling in YBa2Cu30g+x and related materials, and 
use the results to interpret measurements of the optical magnon energy in 



YBa2Cu306.2 and of the NMR "cross-relaxation" rate in Y2Ba4CuyOi5. We 
estimate J± ~ 14 meV in both materials, and Xmax/Ms ~ 100 states/eV-Cu in 
Y2Ba4Cu70i5 at T = 100K. We show that there is at present no obvious con- 
tradiction between the results of a widely-used analysis of NMR experiments 
and the results of neutron scattering experiments in the Y-Ba system. 

We argue that the 41 meV excitation observed in superconducting 
YBa2Cu3C>7 is a collective mode pulled down below the superconducting gap 
by interactions, and that the observed antisymmetry under interchange of 
planes follows from the non- negligible value of J±. 
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I. INTRODUCTION 



In the yttrium-barium (Y-Ba) family of high temperature superconductors the basic 
structural unit is a "bilayer" consisting of two CuC>2 planes; the bilayers are separated by 
CuO chains and the coupling between bilayers is very weak. Neutron scattering)! and more 
recently NMR experiments^ have shown that the Cu spins on adjacent planes in a bilayer 
are coupled. Intra-bilayer coupling has been shown theoretically to lead to a "spin gap" 
similar to that observedi in NMR experiments on underdoped members of the Y-Ba family. 
In view of the great importance of the spin gap phenomenon, a quantitative analysis of 
the spin dynamics of a bilayer is desirable. In this communication we provide this analysis 
and use it to interpret NMR and neutron scattering experiments. Our new results include 
(a) a calculation of the quantum renormalization of the optic spin-wave gap in the S = \ 
Heisenberg bilayer, (b) the proper analysis of the NMR "cross-relaxation" experiment!, (c) 
refined estimates of the inter-bilayer coupling and in-plane antiferromagnetic susceptibility, 
(d) a demonstration that there is no obvious contradiction between the results obtained from 
an analysis of the NMR data in terms of the "MMP" mo deffl and the results obtained 
from neutron scattering^ and (e) an improved understanding of a theory for the 41 meV 
peak observed^ in neutron scattering experiments on superconducting YBa2Cu307. A brief 
account of this work has been presented elsewhere^ 

Focus on the two planes of a bilayer, and neglect coupling to other bilayers. Label the 
spin degrees of freedom by an index a = 1, 2 distinguishing planes, and a site index i. 
Consider the susceptibility 

X ab (q,oo) = j°° dt e i"t+H^-*i)) ^(0)]} . (1) 

Here 5* is the usual spin | operator. Because of the symmetry under exchange of planes, 
X ab has only two independent components, x n (g,u;) = x 22 (g,u;) and x 12 (g,u;) = x 21 (g, u;). 
The two independent components of \ may be taken to be the even and odd (under the 
interchange of planes) components Xeven,odd = Xn ± Xi2- The main goals of this paper are 
to relate experimental observations to x 11 an d X 12 an d these to Hamiltonian parameters. 



The rest of this paper is organized as follows. In section II we derive an expression 
for the optic spin wave gap for the bilayer Heisenberg model and use the result to obtain 
from data0S an expression for the bilayer coupling Jj_ in YBa2Cu30g.2- In section III we 
formulate, solve and apply to data a model for the spin dynamics of a metallic (but non- 
superconducting) bilayer, obtaining estimates for J±, the antiferromagnetic susceptibility, 
and the temperature dependence of the correlation length. In section IV we extend the model 
to the superconducting state and show that it accounts naturally for the observation of the 
superconducting state 41 meV peak observed in YBa2Cu3O70 and the lack of measurable 
signal in the normal state. Our results provide a new interpretation of a previously published 
numerical calculation!!!, explain why the effect was not seen in another calculatiori0 and 
show that an objection raised by other workers^! is not valid. Section V is a conclusion. 

II. BILAYER HEISENBERG MODEL 

In this section we estimate the optic magnon energy of the bilayer Heisenberg model 
defined by the Hamiltonian 

H = lj y E$ a) -3& + J± y E$S! (2) 

Z ii i 
a 

Here a = 1,2 labels planes, i labels sites of a two dimensional square lattice, 5 labels the 
four in-plane nearest neighbors of a site, and the Si are the usual S = \ spin operators. 

We solve Eq. |] via the Schwinger boson method!, which has been shownll to give very 
accurate estimates for the quantum renormalization of the single-plane S = ~ Heisenberg 
model. We previously used the Schwinger boson method to analyze the quantum disordered 
phases of Eq. The analysis presented here amounts to an extension of that work to the 
ordered phase; we therefore omit all details of the derivation of the equations, and mention 
only necessary notation and new features. 

In brief, to use the Schwinger boson method one introduces Bose fields so that 
£>ia = b]a^ ] a a pbtf) . The fields are subject to the constraint (25 + 1) = Z) a &la Ka- The 



S ■ S interactions become four-boson interactions, which one decouples with Hubbard- 
Stratonovich fields Q (for the in-plane interaction J) and A (for the between planes in- 
teraction Jjl). One must also introduce a Lagrange multiplier /x to enforce the constraint. 
The mean field equations are given in Eqs. 2.9 of Refill and are rewritten here for conve- 
nience: 

d 2 k n 



(2tt) 2 ^ 



d 2 k (Q 7fc + A) 7fc 
(2vr) 2 oo k 

d 2 k (Q 7fc + A) 7fc 
(2tt) 2 



1 + 26(f)" 


= 1 + 2S 




1 + 2b(^)_ 


= Q/2J 




1 + 


= 2A/J ± 


(3) 



Here 7 ^ = — ^(cosk x + cosk y ),uj k = y/i 2 — (Q 7 fc + A) 2 and b(x) = (e x — In the 

limit J_i_ — > these equations reduce those studied by Arovas and Auerbachll!, except that 
in our conventions Q is four times larger than in theirs. At low T for parameters such that 
the ground state is ordered at T = one has 



fi 2 = (Q + A) 2 + (k/2) 2 
with inverse correlation length k given by 

k = K exp(—27ip s /T). 



(4) 



(5) 



Equation | defines the spin stiffness p s . 

In the T — > limit the integrals in Eqs. |3] may be evaluated up to corrections of order 
T by noting that the terms proportional to b(uj k /T) are dominated by the divergence as 
ujk —> while in the terms proportional to 1 one may set k = 0. If T is less than the optic 
mode gap one finds 

Q = 2J[2S +l-I oa + I ob ] 
1 



A 



J ± [2S 



Ps = -Q[2S + l-I 



(6) 



Here I oa ,b,c are respectively the T = values of the first, second and third integrals on 
the left hand sides of Eqs. |[ They depend upon the ratio A/Q and in the limit A <C Q are 

2 [A fA\ 

/c=i.39--^ +0 y 

2 [A {A\ 

Expressions for dynamic susceptibilities are given in Appendix C of Ref.0. In the limit 
of interest here the absorptive parts are dominated by poles representing undamped spin 
waves. In the acoustic sector the dispersion is given by u = u>k with k near Q = (ir, ir). In 
the optic sector the dispersion for wavevectors near Q is given by uj = uik+Q- In particular 
the acoustic spin wave velocity is c 2 = \Q{Q + A) and the optic-mode gap u opt = 2\/AQ. 
Neglecting terms of relative order A/Q we find c = 2\^2Z C (S)JS and p s = 2Z Pa (S)JS 2 with 
Z c (jj) = 1.16 and Z Ps {^) = 0.71, in agreement with previous worktilB. [The factor of two in 
p s occurs because each unit cell contains two Cu ions, one in each plane; for T oj opt the 
two planes move together]. 

The equation for A may be rearranged to read 

j± = — — ( 8 ) 

If one takes the limit J±/ p s — > and sets Z Ps = 1, this reduces to the usual uj opt = 



A^JJlS. 

Recent neutron scattering measurements^'^ have reported u op t = 70meV for 
YBa 2 Cu306.2- Use of the canonical values J = 100 meV and Z Pa = 0.7 implies J± = 14meV. 
Our value for J± differs from that quoted in the experimental papers because those works 
assumed Z Ps = 1 and neglected the second term in the denominator of Eq. Some caution 
is necessary in using this value because the measurements were made on samples of nominal 
composition YBa 2 Cu306.2- The effects of the extra 0.2 of oxygen might lead to substantial 
changes of Z Pa from the ideal-system Z Pa ^0.7 value. 



We can also obtain from the formalism an expression for the strength of the spin-wave 
poles. In the regime of interest the momentum integrals are dominated by the infrared 
divergence at Uk <C T; we find, for the susceptibility per plane, 

x"(q, u) = Wf>s [8(u - uj q ) - 8(u + u q )\ (9) 

UJ q 

For the acoustic branch, u q = cq; for the optic branch, uj q = uj opt . The prefactor may be 
written as AnSV2Z Ps /Z c (qa) or (using c 2 = p s /x±) 4>irSy/2Z x Z c / \qa), where a is the lattice 
constant. This susceptibility is, up to an apparently omitted factor of Z C) consistent with 
that given in Eq. 8 o Those authors define x as j-(S S + ) and measure it per bilayer. 
The factor of two from S~S + versus S Z S Z and the factor of two from bilayers leads to a 
difference of ir/2 in our conventions. 



III. METALLIC BILAYER 

This situation is more complicated and less well defined than the insulating antiferromag- 
netic case because there is neither a generally accepted theory of the magnetic dynamics of 
a single layer nor a generally accepted theory of the interlayer coupling. One must therefore 
proceed phenomenologically. Despite the obvious limitations, we believe that this is worth 
doing because the spin dynamics of metallic cuprates are a subject of continuing interest 
and, as we shall show, recent measurements of interplane effects yield information about 
both the between-planes coupling and the in-plane spin dynamics. 

Further analysis requires a model. We shall assume that the magnetic dynamics of 
a single plane may be described by a susceptibility Xoilj^) whose form we discuss further 
below. We also assume that the only interplane coupling is the magnetic one, J±y^sf^ -Sf . 

i 

This is not an important restriction. Retaining the off-diagonal term Xq 2 m our formalism 
leads only to a modest renormalization of J±, as shown also Finally, we assume that 
the effects of J± may be modeled via the RPA. Thus we write 
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x- 1 



-J± Xo 1 ^^) 



(10) 



We further assume (as has been done in previous analyses of NMR in high T c material: 
that Xo(q,u) has the scaling form Xo(q, = Xo£ 2 f /£ z ), where / is a scaling function 
normalized so that /(0, 0) = 1, z is the dynamical exponent, £ is a correlation length 
measured in units of the lattice constant and q is measured from an ordering wavevector 
Q which for the present discussion is arbitrary. Q is believed to be of the order of (tt, tt) 
in high T c materials. From Eq. [L0] we see that x n (g,a;) = Xo{q.i u )/0- ~~ (J±Xo(Qi uj )) 2 ) 
and x 12 (?) w ) — J±Xo(q, w ) 2 /(l ~~ (J±Xo(q, UJ )) 2 )- The crucial parameter controlling the 
susceptibilities in the static limit is 

I± = J±Xmax = J±X0^ 2 - (11) 

Here Xmax is the maximum value of the in-plane spin susceptibility. We must assume I± < 1 
so that the material has no long range order. If I\ <C 1 then x 11 ~ Xo and x 12 ~ J±Xo- l n 
this limit the interplane coupling has a weak effect and the RPA (which in this limit is just 
a perturbation expansion in J±) is an appropriate model. On the other hand, if I\ « 1 then 
the interplane coupling is strong and the use of the RPA may be questioned. 

We now turn to the NMR experiments of interest. These are T2 experiments performed 
on Y2Ba4CuyOi5, a material in which the single-chain structure of YB^CusCV alternates 
with the double-chain structure of YP^CuiOgBS. As a result, atoms on different planes of 
a bilayer have somewhat different local environments and therefore somewhat different NMR 
resonance frequencies, which may be independently studied. Despite the differences in local 
environment, the electronic properties of the two planes seem essentially identical - in par- 
ticular, ratios of Knight shifts and in-plane relaxation rates are temperature-independent@. 
We therefore believe that the differences in observed relaxation rates and Knight shifts are 
due to differences in hyperfme couplings. An alternative view is that the electronic suscep- 
tibilities on the two planes differ. This difference may be included in our formalism and 
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analysis, and will not change the results in an important way (basically, one replaces x 11 by 

(x iy 2) i /2) _ 

Now the NMR T 2 measures the rate at which a nuclear spin is depolarized by interacting 
with other nuclear spins, i.e. it measures the nuclear-spin-nuclear-spin interaction strength. 
In high T c materials the dominant contribution to the nuclear-spin nuclear-spin interaction 
comes from polarization of electronic spins, and may be related to the static limit of the real 
part of the electronic spin susceptibilityil. In Y 2 Ba 4 Cu70i5 it is possible to measure T 2 , the 
rate at which a spin in one plane is depolarized by spins in the same plane, and T 2 j_, the 
rate at which a spin in one plane is depolarized by spins in the other plane. T 2 is related to 
the electronic spin susceptibility by@ 



T. 



(«) 



EW 0) V] 2 - E4 a) V 



1/2 



(12) 



while T 2 _|_ is given byll 



1 



T, 



2± 



E KM 



2) 12 
q A 



1/2 



(13) 



Here we have allowed for the different hyperfine couplings in the two planes observed 
experimentally. We have calculated T 2 and T 2 j_ from Eqs. ([H],|12|,[13|) . The precise values 
obtained depend upon the form chosen for /((?£, u — 0). We have used two forms for f(x) = 
f{q^,uj = 0): a Lorentzian, f(x) = 1/(1 + x 2 ), and a Gaussian, f(x) = exp(— log(2)x 2 ) (the 
log(2) is introduced so f(x = 1) = 1/2). We measure ^ in units of the lattice constant, we 
set h = 1 and assume that the hyperfine couplings can be approximated by their values at 
Q, Aq. We find 

1 



T, 



(14) 



and 



A { $A%J ±x ieg^L 



T, 



2_L 



(15) 



where <?j n and g±_ are defined in terms of the function f(x) via 



1 

2^ 



dx - _„ — 7-t ~ I / dx 



o - {I -IIP) 2 2tt£ 2 \Jo I -IIP 



(16) 



1 f°° f 

^ ) = 27io dx {i-iipy (17) 



In writing Eqs. [T6l,|TT| we have assumed that the correlation length is so long that lattice 



effects may be neglected. We have investigated this issue by performing the exact integrals 
numerically The parameter governing the size of the lattice effects is (7r£) -1 ; for ( > 1 we 
have found that lattice effects are negligible. 



In Fig. 1 we present the calculated results for sJT^T^ /T 2 ± = I±g±{I±)/gi n {I±)- From 



the experimental values! \/t 2 (1) T 2 (2) /T 2± = 0.15 at T = 200 K and \/t 2 (1) T 2 {2) /T 2± = 0.30 at 
T = 120 K we obtain I ± « 0.2 at T = 200 K and I ± w 0.4 at T = 120 K. Note that even at 
the lowest temperature we find that the system is in the small I± regime in which T 2 /T 2 ± is 
linear in I±, suggesting that the between-planes coupling is sufficiently small that the RPA 
formula is justified. 

The magnitude of J± may be determined if Xmax is known and conversely. In the NMR 
and bulk susceptibility literature x is defined in terms of the \gfMsS, g^BS] correlator, and 
results are presented as \l A*b- This quantity differs by a factor of four from the susceptibility 
defined in this paper; we refer to it as -^ ex P. If the susceptibility were only weakly g-dependent 
then the measured uniform susceptibility Xmfform ~ 2 states/eV-CuS would provide a good 
estimate for 4x max and our value I± ~ 0.4 would imply Jj_ w 0.8 eV. Such a value is very 
difficult to justify on microscopic grounds because the insulating antiferromagnetic parent 
compounds of the high T c superconductors have in-plane exchange constants Ji n _ p i ano ~ 0.12 
eV and, as shown in the previous section, Jj_/ Jin-pianc ~ 0.1. It seems very unlikely that 
doping would lead to a Ji > ^in-piane- Therefore, we believe the cross-relaxation results 
imply Xmax 3> Xuniform m metallic, superconducting materials. A similar conclusion has 
been drawn from analyses of in- plane NMrISS, but these analyses required assumptions 
about magnitudes of hyperfine couplings and correlation lengths. These quantities drop out 
of the present analysis. The estimate J± ~ 10 — 20 meV has been obtained from band 



structure calculations!!, implying Xmax ~ 80-160 states/eV-Cu, while in Ref.El the estimate 
J_i_ « 25meV was presented. 

We now turn to a more detailed discussion of the value of Xmax inferred from the dis- 
cussion of the in-plane NMR. This discussion is based on the results contained in the recent 
paper of Barzykin and Pine which summarizes and analyzes a wide range of NMR data. 
(The results proposed in Ref.@ lead to corrections which are smaller than our estimated 
uncertainties so we have not incorporated them). This paper is extremely useful, but has 
a few limitations: in particular, explicit error estimates are not given, and the analysis for 
YBa 2 Cu306.63 does not include interplane effects argued below to be important. We believe 
the uncertainties in values for YBa 2 Cu30 6 . 63 are substantial and difficult to estimate with 
accuracy. 

Barzykin and Pines show that once the hyperfine coupling is determined, the measured 
in-plane T 2 _1 rate is proportional to the quantity Xmax/C Their Table I and Fig. 2 gives 
for YBa 2 Cu 3 7 the estimate Xq P {T = 100K)/£ =30 states/eV-Cu. We believe 10% uncer- 
tainties in hyperfine couplings are reasonableil, implying 20% uncertainties in Xmax/£- To 
obtain an estimate for Xmax & n estimate for £ is required. We arguedil that the uncertain- 
ties in the NMR estimate for £ are large; we suggested, on the basis of earlier data, that 
for YBa 2 Cu3C>7 at T=100K a plausible estimate is 0.5 ~ £ ~ 2 with the most likely value 
~ 1 — 1.5@. From Fig. 14 and Fig. 16 of Ref.0, which are based on more recent data, one 
obtains the larger estimates £/a = 1.9 — 2.3. For YBa 2 Cu30g.6 it seems clear that £ must be 
larger than in YBa 2 Cu3C>7; The analysis of Ref.0 yields £/a(T = 100K) = 7.4. This value 
seems to us too large, but regardless of the precise value of £, it is clear that the larger T^ 1 
and presumably larger £ implies that Xmax in YBa 2 Cu30g.63 is larger by a factor at least 
of order four than in YBa 2 Cu 3 7 We thus roughly estimate that Xmax i n YBa 2 Cu 3 7 ~ 50 
states/eV— Cu (with ~ factor of two uncertainties) and Xmax i n YBa 2 Cu306.63 is at least 
250 states/eV-Cu (with rather larger uncertainties). The value obtained from the results 
in Ref.0 is Xmax ~ 440 states/eV-Cu. As previously remarked, Y 2 Ba4Cu70is has a doping 
intermediate between the two materials, but probably closer to O7. The previously derived 
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estimate Xmax ~ 80 — 160 states/eV-Cu is thus consistent with the results of the in-plane 
NMR, so it is likely that the effective between-planes coupling does not change much with 
doping. 

The rapid dependence of Xmax with doping inferred from in-plane NMR measurements 
suggest that for YBa 2 Cu30 6 . 63 and lower dopings, Xmax is so large that I± ~ 1. For these 
materials the planes are thus very strongly coupled at low energies. The simple RPA formula 
we have used is unlikely to be accurate, and the between-planes coupling must be taken into 
account in interpreting even in-plane NMR experiments. A preliminary attempt along these 
lines has appeared^. 

In summary, the cross-relaxation experiment shows that the real part of the suscepti- 
bility at some non-zero wavevector q is much larger than the uniform susceptibility. The 
temperature dependence of the T 2 rates must be due to the temperature dependence of 
this antiferromagnetic maximum. Two scenarios have been proposed for the temperature 
dependence: in the antiferromagnetic scenario the temperature dependent quantity is the 
correlation length £. In the generalized marginal Fermi liquid scenario the temperature 



regime where T 2 ± is linear in I± the antiferromagnetic scenario predicts T 2 3 /T 2 ± is tempera- 
ture independent, while the marginal Fermi liquid scenario predicts T 2 2 /T2± is temperature 
independent. The experimentally determined ratios are plotted in Fig. 2 and are more 
consistent with the antiferromagnetic scenario. 

The imaginary parts of the two independent susceptibilities Xeven an d XoM are measur- 
able via neutron scattering because they have different dependences on q z , the momentum 
transverse to the CuC>2 planesQ. Neutron scattering experiments have been performed on a 
variety of metallic members of the yttrium-barium family of high-T c materialsS9Iil. The 
experimental result is that only Xodd is seen. At frequencies less than 30 meV and tempera- 
tures less than room temperature the even parity fluctuations are claimed to be completely 
frozen out. The theory of neutron scattering in high T c materials is presently controversial. 
There is no generally accepted model which correctly accounts for the observed lineshapes 



dependent quantity is the overall amplitude 
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and temperature dependences. To investigate the connection between the cross-relaxation 
experiments and neutron scattering we have chosen to calculate the ratio of the g-integrated 
even and odd parity susceptibilities. This ratio is insensitive to the precise details of the 
susceptibilities. For defmiteness we used the "MMP", dynamical exponent z = 2 ansatz 
Xo(<??^) = x/(£~ 2 + q 2 — ioj/T). Here T is a microscopic spin relaxation time. The results 
depend on I± and on ujsf = r/£ 2 , which we define here to be the softest spin fluctuation 
frequency of a single plane. Of course the between-planes coupling will reduce this frequency 
for the odd parity channel and increase it for the even channel. Results are shown in Fig. 
3 for several values of I±. We see that the relative weight of the even parity fluctuations 
becomes small only for I± > 0.5. We believe that the neutron results, which seem to require 
a I± > 0.5, are not in contradiction to our analysis of the cross-relaxation experiment, which 
yielded a I± < 0.4, because the strongest neutron evidence for locked bilayers was obtained 
from a study of YBa2Cu3O6.50, which as we have previously noted is much closer to the 
magnetic instability than Y 2 Ba4Cu70i5, and therefore may be expected to have I± ~ 1. 

As a side remark we note that a recent neutron scattering work provided a bound on 
the magnetic scattering in YBa2Cu3O70. Bearing in mind the factor of (g 2 ir/2) difference 
between NMR and neutron conventions discussed above, we note that the upper bound 
implied by the neutron experiment for Y Ba^Cu^Oj is Xmax < 300 states/eV-Cu. This 
bound is comfortably larger than the NMR estimates obtained by interpreting in-plane 
NMR experiments. At this time, therefore, there is no evidence of a contradiction between 
neutron scattering experiments and the "MMP" scenario for interpreting the NMR. 



IV. SUPERCONDUCTING STATE 

Another experiment relevant to bilayer coupling in the Y-Ba system is the observation of 
a sharp peak at an energy of 41 meV in the antisymmetric susceptibility of superconducting 
YBaaCuaOji. A corresponding peak in the symmetric channel has not been observed00, 
and an upper bound on the intensity in the symmetric channel of about 30% of the inten- 
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sity in the antisymmetric channel has been established^. A number of theoretical works 
have investigated this peak (see, e.g. roil and references therein). We believe that the 
most likely explanation is that the peak is a collective mode pulled down below the su- 
perconducting gap edge by interactions, and that its appearance only in the antisymmetric 
channel is due to the non-zero J±. This explanation is essentially that given by Liu et al.0 
(although they did not interpret their results - which were basically numerical - in precisely 
this manner). An objection to this interpretation was raised by Mazin and Yakovenko0, 
who argued that within such theories an observable peak should also exist in the symmetric 
channel. In the remainder of this section we present a simple RPA theory of the 41 meV 
mode. Although the RPA is presumably not quantitatively accurate it is in fact the basis 
of most previous work (@-@) and is to some extent analytically tractable. Our treatment 
gives insight into the results of0, shows why the effect was not seen by0 and shows that 
the objection oil is unfounded. 

The fundamental object in the RPA is the bare electron polarizability Xoo- This is given, 
at T = in a superconductor, by e.g. Eq. 1 of Refill (although one must multiply their 
expression by 2 to convert to our conventions). In two spatial dimensions the constraints that 
the initial and final states lie on the Fermi surface completely specifies the kinematics (up 
to discrete lattice symmetry operations). Thus the threshold behavior of x" f° r frequencies 
near the gap edge may be found analytically. For Q — (tt, tt) and bearing in mind that 2Q 
is a reciprocal lattice vector one finds 



(n \ reg . 4 \l |A po A po+ 



~ v Po v Po+Q sin 6 



lot 



uj - A 



A 



Po- 



reg 

Xoo 



1 + A log 



A PO A PO+C 



u-A P0 - A P0 . 



(18) 



Here p is a wavevector such that e po = e P0+ Q = 0, 9 is the angle between v Po and 
v Po+ q and it is assumed that the Fermi surface does not pass through a van Hove point 
(so v po 7^ 0) and 6 ^ 0, tt. Xoo is a function which is not singular at the gap edge; one 
expects Xoo to be basically equal to the normal state susceptibility, and A ~ (A/Ep). Of 
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course, as the Fermi surface approaches the van Hove point, A increases. The dispersion 
e p = —2t(cosp x + cospy) +4t' cosp x cosp y — fi with if — t/2 and \x = — 1.5t has been claimed 
to describe YE^CuaOyB For this dispersion Xoo — 0.6/t and A = 0.8(A/t). Thus using 
t ~ 0.2eV and A ~ 20 meV, we estimate A ~ 0.1. The coefficient of the logarithm is quite 
sensitive to details of band structure, thus the significance of this estimate is that it may be 
reasonable to assume A is small, but not absurdly so. Note, though, that a key feature of 
the calculation is that at po and po + Q the superconducting gap is at or near its maximum, 
so that A is maximal and there is no quasiparticle damping. Other workers (includingEl) 
have studied a model with t' — 0. For this model, po and po + Q are near gap minima, so 
A is very small and the damping is large, so the effects we will discuss below will not be 
observable. Thus it may be concluded that within the present theory the 41 meV peak is 
very sensitive to details of band structure and gap anisotropy. 

The RPA formula ^( a ' s ) = xoo/(l — {J ± J±)xoo) an d Eq. [18] thus imply that the 
susceptibility diverges at frequencies u a,s given by 



2A P0 - u a > s = 2A P0 exp 



(19) 



(J ± J±)XmaxA_ 

Here Xmax = Xoo/ (1 — ^Xoo) and I± = J±Xmax have the same meaning as in section 3. One 
sees immediately that if the system is not too far from a magnetic instability the large value 
of JXmax m ay compensate for the small value of A. Also, the argument of the exponential 
differs for the symmetric (1 + I±)/(J — J±) and antisymmetric (1 — I±)/(J + Jj_) cases. 
If JXmaxA < 1 and I± is not too small, then only the antisymmetric pole is appreciably 
removed from the gap edge. Finally, the weight in the resonance, J* = / dux , may be 
written approximately 

/a , s = 4tt(2A P0 - u a<8 ) ^ 
J A 

The weight scales with the difference in energy from the gap edge, and is thus much 
smaller for the symmetric mode. It is enhanced by the factor 1/A ~ 10. Because precise 
values for the various quantities are not available it is difficult to compare our estimate 
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of the weight precisely to experiment. In Ref.Ej the value J ~ 80 meV was useda. The 
conventions are those of this paper, so to convert I to the units used inO one must multiply 
by 2/tt. The result is I a,s ~ 40 (1 — uj a,s ) 2A Po . It thus seems that weights of order unity or 
larger in the antisymmetric channel are not difficult to obtain. 



V. CONCLUSIONS 

We have given expressions which may be useful in interpreting experiments relating to 
bilayer coupling in YBa 2 Cu30 6+x . We have determined the quantum correction to the optic 
magnon gap in the insulator, and given the theory of the NMR cross relaxation T 2 exper- 
iment. We have shown that a J± ~ 14 meV is roughly consistent with all experiments, 
and that the cross-relaxation experiment implies that x has a very substantial antiferro- 
magnetic peak, which however is not large enough to have been observed in recent neutron 
experiments. 

We have shown that a J± ~ 14 meV is consistent with neutron scattering experiments 
and with an extension to the bilayer system of the widely - used "MMP" analysisi'0 of 
NMR. We have also noted that there is at present no obvious contradiction between neutron 
scattering and the MMP scenario for NMR. We have proposed an explanation for the 41 meV 
peak observed in the superconducting state of YB^CusOy and have argued that in the Y-Ba 
system at dopings corresponding to YBa2Cus06.63 interplane effects will have a significant 
effect on the low-energy magnetic dynamics of the underdoped materials. Understanding 
the effects is a crucial problem for future research. 
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FIGURES 




FIG. 1. Ratio of cross-relaxation rate 1/T 2 ± to in-plane relaxation rate I/T2 plotted versus 
coupling parameter I± = J±Xmax for Lorentzian (dotted line) and Gaussian (solid line) form factors 
and calculated from Eqs. (6-9). 
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FIG. 2. Experimentally determined ratio of 1/T2j_ to nth power of I/T2 for n=l (□), 2 (•), 3 
(o) in arbitrary units. That the n=3 (o) curve has less temperature dependence than the n=2 (•) 
curve suggests the existence of a growing magnetic correlation length. 
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FIG. 3. Calculated ratio of (/-integrated odd-parity neutron absorption to q-integrated even 
parity neutron absorption, plotted versus frequency for 7j_ — J iSKmax — 

0.0, 0.1 ...0.9. I± = 0.0 

corresponds to the lowest curve and I± = 0.9 to the top curve. 
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